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Abstract 

The possibility of QQqq heavy-light four-quark bound states has been analyzed by means of the chiral 
SU(3) quark model, where Q is the heavy quark (c or b) and q is the light quark (u, d or s). We obtain a bound 
state for the bbhh configuration with quantum number = = and for the ccrih. (J^ = = 0) 
configuration which is not bound but slightly above the D*D* threshold (n is m or (i quark). Meanwhile, we 
also conclude that a weakly bound state in bbnh system can also be found without considering the chiral 
quark interactions between the two light quarks, yet its binding energy is weaker than that with the chiral 
quark interactions. 
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I. INTRODUCTION 



A great interest in studying four-quark (4q) states has been triggered by the observations of 
some very narrow resonances, e.g., D*^(2317) by BABAR ilD, D,ji2460) by CLEO [2], X(3872) 
by BELLE |3], and 7(4260) by BABAR J^]. Because almost all of these new states have the 
same special properties, which are narrow and their structures are difficult to be explained as the 
simple quark-antiquark pair, much attention has been paid and various theoretical explanations 
were proposed. Mostly, these new states were suggested as good candidates of 4q states [5]. In 
this work we will try to see if there are possibilities for the existence of the four-quark states with 
two heavy quarks which is also a significant problem in the new hadron state studies. 

The possible existence of 4q states composed of two heavy quarks (anti-quarks) and two light 
anti-quarks (quarks) has been suggested and studied in the past few years in difl'erent ways |^ 7, sj]. 
For example, the work of Carlson et al. based on the MIT bag model concluded that the state 
bbiid with = 1^ is bound and the ccud is not bound. A different conclusion is obtained by 
Pepin et al. [SD using a pseudoscalar meson exchange interaction instead of the chromomagnetic 
potential, and their results indicated that such interaction binds the heavy-light 4q systems both 
forQ = b and Q = c. Therefore, whether the bound state exists in heavy-light 4q systems is still 
an open question. 

In this work, we focus our attention on the heavy-light 4q states and study their masses in 
the framework of chiral quark model (CQM). Since the heavy-light four-quark state contains two 
light quarks (anti-quarks) which are regarded as constituent quarks, the non-perturbative QCD ef- 
fect should be considered between them. On the other hand we notice that the chiral quark model 

[l^ based on the idea that between the chiral symmetry breaking scale and the confinement 
scale, QCD may be formulated for the light-quark sector as an effective theory of constituent 
quarks interacting through Goldstone boson modes associated to the spontaneous breaking of chi- 
ral symmetry, and CQM has obtained great success in reproducing the energies of the baryon 
ground states, the binding e nerg y of deuteron, the nucleon-nucleon phase shifts, and the nucleon- 



hyperon cross sections jld 



111] . Thus one can expect that the CQM is an effective approach to 



study the 4q systems with light quarks. In this work, the four-quark states with quark content 
cess, cchh, bbss and bbhh are studied in the chiral SU(3) quark (CSQ) model. 

The parameters for the light quark pairs are taken from our previous work, and for the 
heavy-light pairs they are obtained from fitting the masses of (D, D*,Ds, D*, rjc, J/i/r, hdlp)), and 
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(B, B*, Bs, 5*, T]h, T(l5)), the details will be introduced in the section n.B. By using these parame- 
ters, we obtain a bound state for bbhh configuration and no bound state for ccqq. 

Furthermore, we calculate the masses of QQhh without considering chiral quark interaction 
(CQI), and compare them with the masses of QQhh considering CQI. Through comparison we 
uncover that a bound state in bbhh system can also be found without considering CQI, yet this 
bound state is not as strong as the bound state with CQI. 

The paper is arranged as follows. In the next section, the theoretical framework of CSQ model, 
the determination of parameters and the wave functions for QQqq stats are briefly introduced. The 
calculated results are listed in Section in, where some discussions are made as well. Finally, the 
summary is given in Section IV. 

II. THEORETICAL FRAMEWORK 
A. The model 

In the chiral SU(3) quark (CSQ) model, the Hamiltonian of a QQqq four-quark system can be 
written as 



where m, is the mass of the ith quark. To is the kinetic energy operator for the center-of-mass 
motion. Since the light quark is treated as constituent quark, and the non-perturbative QCD ef- 
fect between them has to be considered, thus for the light quark pairs, besides the confinement 
potential and OGE interaction, the chiral quark interactions should be included also, but which is 
unnecessary to be considered for the heavy-heavy and heavy-light quark pairs. The expressions of 
the interactions are written as follows: 




(1) 



,-=l,2,y=3,4 




(34), 



(2) 



For the heavy quark pair {QQ), the confinement potential Vqq (12), which provides the non- 
perturbative QCD effect in long distance, is taken as linear form 



(3) 
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Vqq^ (12) is the one-gluon-exchange interaction, whose expression is 

V.T(12) = 7^.fta; • 4) {- - ^( A * T + i^i^Wr,.) I . (4) 

4 [rn 2 m\ 3mim2 J 

For the light anti-quark pair {qq), we can obtain the forms of V^^*^^ and V^^"^ by replacing the 
/I; • /I'^i in Eq. (3) and Eq. (4) by A'^* ■ /l", and the chiral quark interaction between two anti-quarks 
V|| has the same form (only ■ A^ is replaced by A^* • A^*) as that of two light quarks V^^. In the 
CSQ model, V^^ includes the scalar boson exchanges and the pseudoscalar boson exchanges, thus 

8 8 

V2(34) = J] K„(r34) + 2 V,,SrM), (5) 
and the expressions of these potentials are 

VsMsa) = -C(gch, A)(Af ■ AT)X,(m,^, A, r,,), (6) 

2 

m 

VpsSr34) = C(gch, mps^, A)(^r " -^4*) Y^;^^2K.<„ A, r^,){a^ ■ a,), (7) 

where 

C(g,,,m,A) = ^-^, (8) 
An A^ - m'^ 

Xi(m, A, r) = y(mr) - -y(Ar), (9) 

m 

X2(m, A, r) = y(mr) - {-fY^Kr), (10) 

m 

F(x) = V^ (11) 

and ms^inipsj is the mass of the scalar (pseudoscalar) meson. 

The interactions of quark- antiquark pair (Qq) include two parts: direct interaction and annihi- 
lation part 

y^r ^ yconf ^ yOOE 



The contribution of annihilation part is neglected in this work, and the expression of Vq^ can be 
^fg. Concerning V^-"^ and Vg?' 



obtained from Vq^. Concerning Vq^-^ and Vq9^, the transformation from Vg^ (the expressions of 



Vq^^ and Vq^^ are the same as Eqs.(3) and (4)) to VQg is given by A': ■ A'j — > -A? • A^, while for 



V^J;^ the alteration is A" ■ A" ^ A" ■ A"*. 



B. Determination of the parameters 



The interaction parameters include the OGE coupling constant gt, the confinement strengths 
(aij, alj), and the chiral coupling constant gch- For the light quark pairs, the parameters are taken 
from our previous works 11121 11311 . which gave a satisfactory description for the energies of the 
light baryon ground states, the nucleon-nucleon scattering phase shifts the binding energy of the 
deuteron, where the chiral coupling constant g^h is fixed by 

,2 



4n 25 An 



(14) 



with g^j^ /An = 13.67 taken as the experimental value. The parameters for cq or cc quark pairs 
are taken from Ref. [12], which fitted the masses of D, D* ,Ds, D*^, rjc, J 1^ and hi{\ p). Followed 
the same method, the model parameters for bq and bb quark pairs can be fixed by the masses 
of B,B*,Bs,B*,rih and T(li'). The complete set of the parameters is shown in Table I, and the 
corresponding theoretical results for the masses of Qq and QQ mesons are shown in Table II, 
from which we can see that the theoretical masses of mesons are reasonably consistent with their 
experimental values, respectively. 



TABLE I: Model parameters for the systems with heavy quarks. The masses of u (d) and s quark: m„ ^ 
3\3MeV, m, = AlOMeV. The OGE coupling constant g„ = 0.^ 



nic (MeV) 


1430 


rrib (MeV) 


4717 


gc 


0.58 


gb 


0.52 


ttcu (MeV/fm) 


275 


abu (MeV/fm) 


275 


4u (MeV) 


-155.9 


al (MeV) 


-141.1 


Ucs (MeV/fm) 


275 


Ghs (MeV/fm) 


275 


4s (MeV) 


-124.7 


al (MeV) 


-112 


Qcc (MeV/fm) 


275 


ahb (MeV/fm) 


275 


(MeV) 


-77.8 


al (MeV) 


-39.5 
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TABLE II: The masses (MeV) of Qq and QQ mesons. Experimental data are taken from PDG. 



Mesons 


D 


D* 




D* 


'7c 




h,{\p) 


Exp. 


1867.7 


2008.9 


1968.5 


2112.4 


2979.6 


3096.9 


3526.2 


Theor. 


1888 


2009 


1969 


2130 


2990 


3098 


3568 


Mesons 


B 


B* 


Bs 


B* 


Vb 


T(l^) 




Exp. 


5279.2 


5325 


5369.6 


5416.6 


9300 


9460.3 




Theor. 


5288 


5320 


5317 


5412 


9404 


9460 





C. Four-quark wave function 

For the description of a 4q state (QQqq), the wave function can be written as the following 
form 

^(4^) = MOs')[(QQ)l\cSqq)tcAc' (15) 

where S represents spin, / represents isospin, and C represents color. ifr4y(0s'^) is the orbital part 
and all the four quarks are in (Os)-wave state; [(QQ)'^^ cS^^^'si C2^s c flavor-spin-color part. 



Since the wave function of (QGY^^ a ^^'^ iQ^^Si C2 P^^^^ should be antisymmetric [|l4ll . the possible 



configurations of QQqq are read as follows, 

= (A4,(O/)[((20o,6^(nn)°gJ? (16) 



= (A4,(O/)[(e0,w(«n);,5Jo,io (17) 



= 'A4,(O/)[(e0i,3.(nn)[|3.]° (18) 

^z), = >p4,(0s')[{QQ)am)U^]\i^,(J = 0, 1,2). (19) 

^£ = 'A4,(0/)[(ee)o.6K^^)o,6^]o,i^, (20) 

^F. = 'A4,(O/)[(e0i,3.(^^)i,3^]y,i^,(/ = 0, 1,2). (21) 
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Thus for different quantum numbers (/^; /), the configurations of QQqq states can be written as 



(22) 



r;l) = ^Bfl«J^z)o, 



(23) 



|i";0) = 



and 



(24) 



(25) 



|2";0) = »Ff„ 



(26) 



1) = ^D,. 



(27) 



We solve the Schrodinger equation of the 4q system by using the variation method. The trail 
wave function is taken as an expression of the 4q states with several diff"erent harmonic oscillator 
frequencies w„ 



(28) 



where = —. Then the energies of these states can be obtained. For the 4q states with the same 

1 maii o ^ 



(7^; /), the configuration mixture has been considered in our calculation. 



m. RESULTS AND DISCUSSIONS 

4q state will be stable under the strong interaction if their total energy lies below all the thresh- 
olds of the possible and allowed two-meson decay channels. Therefore, we use the quantity 

AE = EiQQqq) - E„,(Qq) - E„,(Qq) (29) 

to discriminate the stable 4q states. 

We calculated the energies of the 4q states with the quantum numbers: (J^; /) = (0"^; 0), (0"^; 1), 
(l-'lO), (1+; 1), (2+;0), and(2+; 1). The results of ccqq and bbqq are listed in table IE and IV, 
respectively. Our calculation predicts one bound state for the bbhh system with (J^; /) = (1"^; 0) 
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TABLE III: E{ccqq) and AS energies in MeV 



{J' -J) 


(0+;0) 


(0+; 1) 


(l+;0) 


(l+;0) 




(2+;0) 


(2+; 1) 




cess 


ccnn 


ccnn 


cess 


ccnn 


cess 


ccnn 


E 


4402 


4161 


4068 


4462 


4196 


4504 


4231 


^E 


142 


143 


50 


363 


299 


244 


213 


TABLE IV: E{bbqq) and ^E energies in MeV 


{J' -J) 


(0+;0) 


(0+; 1) 


(l+;0) 


(l+;0) 


(i^;i) 


(2+;0) 


(2+; 1) 




bbss 


bbnn 


bbnn 


bbss 


bbnn 


bbss 


bbnn 


E 


10934 


10695 


10576 


10945 


10701 


10959 


10712 


^E 


192 


119 


-32 


216 


93 


135 


72 



and no bound state for the ccqq system, both of which agree with the conclusion that for a ratio 
M/m > 15 (M is the mass of the heavy quark corb,m is the mass of light quark n or s), a collective 
bound state appears, stable against spontaneous dissociation [15]. Meanwhile, we notice that the 
states ccnn and bbnn with (7^; /) = (1+; 0) have the lowest energy than other QQqq states. For the 
(7^; /) = (1^; 0) state, we calculated the energies of the configurations and (Eqs. (16) and 
(18)) and also the configuration mixing between and • The results show that the energy of 
the configuration ¥c is lower and the mixture effect between and *Fc is very small, thus the 
component can be neglected and component \{QQ)\^y{hrifQy]^^ is dominate. In the light 
anti-quark pair {hnf^ y contributes a strong attractive color magnetic force (CMF) which upholds 
Jaffe's di-quark model Jl^, at the same time the pseudoscalar meson exchange interaction of 
{hrif^y is also strongly attractive, thus {hrif^y can be sufficiently attractive for forming a possible 
bound 4q state. On the other hand, in the state \{QQ)\^y{hrifQy]^^ the heavy quark pair {QQ)i^j,c 
contributes a repulsive CMF. When Q = b, the repulsive CMF is much weaker than that of the 
attractive force contributed by the light anti-quark pair component (nn)^^^ because of the large 
heavy quark mass. Yet when Q = c, the repulsive CMF is not weak enough to counteract the 
attractive force. That is why for the same (7^; /) = (1^; 0) state, a bound state in bbqq system can 
be found while there is no bound state in ccqq system. 
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Since CMF and the pseudoscalar meson exchange interaction can contribute enough attractive 
force to form a bound 4q state, it is interesting to find whether the existence of the bound state 
depends on the contribution of CQI or not. Following the above procedure, the parameters without 
considering CQI should be fixed firstly. Here we briefly give the procedure for the parameter 
determination. The one-gluon-exchange coupling constant gu is determined by the mass split 
between and A, and the confinement strengths (a„„,a°„) are fixed by the masses of and A. 
(a„,, a^.J are fixed by the masses of D and D*, while {atu, a^J are fixed by the masses of B and B*. 
Other parameters are taken from previous work iIibIi . The parameters without considering CQI 
are shown in table V, and the theoretical results for the masses of N, A and Qq mesons without 
considering CQI are given in Table VI. 

TABLE V: Model parameters without CQI. The quark masses are as follows: m„^^ = 3l3MeV, nic = 
USOMf'V, mt = AlllMeV. 



gu 0.945 

al, (MeV) -181.68 

a„„(MeV/fm) 187.188 

4^ (MeV) -153.5 

acu (MeV/fm) 275 

al^ (MeV) -138.8 

Uhu (MeV/fm) 275 



By using the variation method in section II.C, the masses of QQqq states can be obtained. 
Here we only show the numerical results for the states cchh and bbhh with {J^;I) = (1"^;0) in 
Table VII. From Table VII, we can see that without considering CQI, the configuration bbhn with 
{J^;I) = (1^;0) is still a bound state and no bound state for cchh configuration. Meanwhile, 

TABLE VI: The masses (MeV) of P, A and Qq mesons. Experimental data are taken from PDG. 





P 


A 


D 


D* 


B 


B* 


Exp. 


939 


1232 


1967.7 


2008.9 


5279.2 


5325 


Theor. 


939 


1237 


1879 


2009 


5287 


5321 



9 



TABLE VII: E{bbhn), E{ccnn) and AE energies in MeV 



Without CQI 


Considering CQI 




ccnn 


bbhh 


ccnn 


bbhh 


E 


4078 


10590 


4068 


10576 


AE 


60 


-18 


50 


-32 



it should be noticed that this bound state is not as strong as that of the case with chiral quark 
interactions. 

In order to see the differences between these two cases, we calculate the contribution of each 
part for the configuration bbhh with (/^; /) = (1"^; 0). We find that without considering CQI, the 
strength of OGE interaction increases and it contributes enough attractive force to form a bound 
state, thus without CQI, bbhh with iJ^;I) = (1^;0) can also be bound. In comparison, when 
considering CQI, the OGE interaction and meson-exchange interaction both contribute attractive 
forces, and the contribution of the pseudoscalar meson exchange interaction is more favored to 
form a bound state. That is why for the configuration bbhh with (/^; /) = (1"^; 0), the bound state 
considering CQI is stronger than without considering CQI. 

IV. SUMMARY 

In this work, we study the masses of QQqq heavy-light 4q system in the CSQ model, and try to 
see the possibilities of the existence of bound states. Our calculation shows the existence of only 
one bound four-quark state with (7^; 7) = (1"^; 0) in the bottom sector and no bound state for the 
charm sector. Besides, by comparing the masses of the same four-quark state while considering 
CQI and without CQI, we find a similar result that bbhh with (J^; 7) = (1"^; 0) is bound but this 
binding energy is relatively weak. 
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